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A g'-LOAD DISPLAY FOR REMOTELY PILOTED VEHICLES AND 


OTHER AIRCRAFT THAT ARE SUBJECT TO MANEUVERING CONSTRAINTS 

James C. Howard 
Ames Research Center 


SUMMARY 


The absence of kinesthetic cues deprives RPV pilots of essential feed- 
back, particularly during high speed maneuvers. A display has been designed 
which presents g'-load information in a sufficiently dramatic form to compen- 
sate for this deficiency. The display embodies all the motion vector compo- 
nents that contribute to maneuvering accelerations. In one mode of operation, 
the device displays the angular velocity components as variable length point- 
ers emanating from the center of a circle whose radius is determined by the 
linear velocity of the vehicle. Inasmuch as this mode of presentation empha- 
sizes the role of the angular velocity vector, while at the same time display- 
ing the influence of the constraint imposed by the linear velocity, it is well 
suited to the application considered in this report. An alternative mode of 
operation permits the pilot to select a display in which g loads are presented 
directly as variable length pointers emanating from the center of a circle 
whose radius is again determined by the linear velocity of the vehicle. Each 
circle defines a safe region of operation, and penetration of the circumfer- 
ence of the appropriate circle by a pointer implies that the structural design 
load has been exceeded . 


INTRODUCTION 


The overall objective of the remotely piloted vehicle (RPV) program is 
to develop an aeronautical research tool that is well suited to aerodynamic 
research; to obtain high angle-of-attack data up to and including post-stall, 
pre-spin conditions, and to assess advanced control systems under these condi- 
tions. In one of its more important applications, the RPV program represents 
an attempt to reduce the cost and the danger involved in testing high- 
performance aircraft. The lack of full-scale flight test data for such air- 
craft is due to the risk involved in testing them. The RPV approach offers 
a very attractive potential for performing large scale flight testing of haz- 
ardous tasks at low risk and relatively low cost (ref. 1). As the RPV program/ 
is presently constituted, remotely piloted vehicles are to be flown from 
fixed-base control centers. In this connection it should be noted that in the 
absence of inertial motion cues, the RPV pilot is more likely than the pilot 
of a manned vehicle to exceed the structural capability of his aircraft during 
high-speed maneuvers. 

It is with this aspect of the RPV program that the present report is con- 
cerned. The g-load display concept proposed is based on a geometrical inter- 



pretation of maneuvering accelerations. Although the maneuvering accelera- 
tions can be depicted in vector form, it is more expedient for present pur- 
poses to reduce the acceleration to scalar form. The most convenient form is 
realized by squaring the acceleration vector and reducing the subsequent quad- 
ratic form to the canonical form by an orthogonal transformation. The advan- 
tage of the canonical representation is that it reveals the geometry of the 
acceleration in a form which is easily interpreted for display purposes. 

It will be seen that gravity can be included without sacrificing the geo- 
metrical simplicity of the maneuvering acceleration display. This is accom- 
plished by expressing the gravity vector in terms of a virtual angular veloc- 
ity vector, which can then be combined with the maneuvering angular velocity 
vector to give the total inertial and gravitational acceleration. 


ANALYSIS 


Maneuvering Accelerations 

Vehicle accelerations resulting from aerodynamic, thrust, and gravita- 
tional forces are given by Coriolis' equation 

= a ( 1 ) 

The geometrical properties of acceleration become more evident when ac- 
celeration is expressed in scalar form. The scalar equivalent of equation (1) 
is 
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Mor^eover, in the absence of accelerations due to d and B, 37/8t is paral- 
lel to 7, and in this case 97/ and co x 7 are perpendicular vectors (see ap- 
pendix). Because of large damping, aircraft response time is short and both 
d and 3 rapidly tend to zero (ref. 2). Hence, the assumption of perpendicu- 
larity is a valid one and gives rise to conservative estimates of accelera- 
tion. This assumption is reflected in the following equations: 
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Now 97/9t is a vector approximately parallel to the flight path, and co x 7 is a 
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vector normal to the flight path. From the point of view of structural inte- 
grity, inertia force components normal to the flight path are more important 
than components parallel to the flight path. This is due to the fact that 
components normal to the flight path subject the wings and empennage to bend- 
ing moments that are reacted by the least moments of resistance of these 
structural units. Hence, the stress-inducing capability of inertia force com- 
ponents normal to the flight path is much greater than that of components 
parallel to the flight path. For this reason the g loads produced by inertia 
forces acting parallel to the flight path will be omitted. It should be 
noted, however, that due to the existence of angles of attack and sideslip, 
the inertia force components normal to the flight path are not always normal 
to the wing and stabilizer plan forms. However, the errors introduced by 
assuming normality are on the conservative side. 

A maneuvering acceleration resulting from changes in aircraft attitude is 
defined as follows: 
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Relative to an orthogonal Cartesian coordinate system, the vector a as- 


sumes the following form: 

- r A A 

= \_(qu) - vv)i + (ru-pw)g + (pv - qu)k] 
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where j and k are a triad of mutually orthogonal unit vectors. 
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As already indicated, it is more expedient for display purposes to reduce 
the acceleration from the vector form of equation (4) to the following scalar 
forms : 


(iiixV) • To) X 7^ = a ^ (6) 

m 

therefore 

a ^ - T7 • (7) 

m 

In terms of the linear and angular velocity components, equation (7) may 
be expanded as follows: 

= {_(v^ +w^)p^ + (u^+w^)q^ + (u^ + v^)v^ - 2(uvpq + uwpr + vwqvJl (8) 

Equation (8) reveals that is a quadratic form. By a suitable trans- 
formation of axes, a quadratic form can be reduced to the canonical form 

aj- = rXiP^ + 

where p, p, r are the components of the angular velocity vector in the new 
system of axes and are the corresponding eigenvalues (ref. 3). 

In the present case, the required canonical form can be obtained by 
choosing a system of axes such that the angular velocity component p coincides 
with the linear velocity vector 7. Although the axes for the specification of 
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the components q and v may be chosen 
arbitrarily, provided they are ortho- 
gonal to each other and to the linear 
velocity vector V, the arrangement 
shown in figure 1 has been adopted. 
Thus, relative to the new reference 
frame, the eigenvalues are 

h = 0 , \2 = , A 3 = 

and assumes the form 

(V^q^ + V^r^) = a 2 (9) 

^ m 


which represents a circular cylinder. 


Figure 1. — Systems of reference axes, 

including body and wind. It follows from equation (9) 

that, for an allowable acceleration 
Uq, the angular velocity vector co is 

bounded by the surface of a circular cylinder of radius CL^IV, having for its 
axis the linear velocity vector V. 


Gravity Components 

In order to retain the geometrical simplicity of the maneuvering g~load 
display, the gravity components may be treated as virtual angular rates. This 
results in adjustments to the pointer length to reflect the instantaneous mag- 
nitudes of the gravity components, while leaving the circle unchanged. For 
this purpose, it is necessary to refer the gravity components to a set of wind 
axes. The axes are chosen such that the x axis coincides with the wind vector, 
and the y and 3 axes assume the orientations shown in figure 1. 


Wind axes are related to body axes as follows: 
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are column vectors of coordinates in the wind and body axes respectively, and 


B = sin ^(v/V) j a=tan ^ (w/u) 

The gravity vector has body axis components, 

_ A A A 

g = g(- sin Qi + cos 0 sin + cos 0 cos cj)fe j 

where 0 and <j) are Eulerian angles relating the body axes to a set of earth 
fixed reference axes (ref. 2). Hence, the components of gravity in the direc- 
tion of the wind axes of interest are: 

g— = g(sin 0 cos a sin B + cos 0 sin 4> cos B - oos 0 cos (|) sin a sin 

and 

g^ = g( sin 0 sin a + cos 0 oos (|) oos a) 

Moreover , 


q = (-p cos a sin Q + q oos Q - r sin a sin 


and 


r = (-p sin o. + r oos a) 

The Q'-load display is required to indicate the instantaneous value of the in- 
ertial and gravitational forces to which the aircraft is exposed. The sum of 
these forces is given by 


therefore 


-ma + mq = - m(a - q) - mf 
m ^ m ^ '' 


f = - (a^-g) = - C^xV-g) 
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A A A 
where i, j , and k 
the directions of 
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are a triad of mutually orthogonal unit vectors specifying 
the wind axis. Therefore, 


= I \g(sin 0 cos a sin B + cos B sin (j> cos B - cos 0 cos <p sin a sin Q) 

- V(-p sin a + r cos aJ] j + [g(sin 0 sin a + oos 0 cos <}i cos a) 

+ V(-p cos a sin Q + q oos ^ - r sin a sin k 


This equation may be rewritten as follows: 


5 



f=v 


^(sin 6 cos a sin g + cos 0 sin (j) cos 6 ~ oos 0 cos ip sin a sin Q) 


- (-p sin a + r cos a) 


0 + 


^(sin 0 sin a + cos 0 cos <j) cos a) 


+ (~p cos a sin ^ + q cos g - r sin a sin gj 


The combined inertial and gravitational accelerations may now be ex- 
pressed in the alternative form 


f = - Q X 7 


where 


_ A A A 
a = Pi + Qj + Rk 


Hence , 


mf = mV 


' A A 
Qk - Hj 


( 10 ) 


therefore 


Q = (- p cos a sin g + g cos g - r> sin a sin g) 


+ sin 0 sin a + cos 0 cos (j> cos a) 


( 11 ) 


or 


Q = (q + q ) 
9 


where 


q = (~p oos a sin ^ + q cos g - r sin a sin g) 


and 


q^ = sin 0 sin a + oos 0 cos cp cos a) 


The term q is the virtual angular velocity relative to the j wind axes cor- 
responding^to the gravitational component of acceleration. 


Likewise, 
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R = 


(- p sin OL + V cos a) - sin 0 cos a sin 0 + cos 0 sin 4> cos 3 


- cos 0 cos (j) sin a sin 


( 12 ) 


or 


where 


R = (r + r ) 
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r = (- p sin a + r cos d) 


and 


r 
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- ^(sin 0 cos a sin 3 + cos 9 sin (f> cos 3 - cos 9 cos <^> sin a sin 


_ A 

where Tg is the virtual angular velocity in the direction of the H axis , cor- 
responding to the gravitational component of acceleration. 


Therefore, the motion vector components to 
be utilized for display purposes are Qj i?, and 
which are related as follows: 

= ^2 ( 13 ) 

The effective angular velocity, which now 
includes the virtual components due to gravity, 
is 


(Q'^ + R^) = = f2/yl ( 14 ) 



_A 

and is inclined at an angle n to the j axis, 
where 


n = tan ^ (R/Q) 


(15) 


Equation (13) has the same form as equation 
(9) . It follows that for an allowable accelera- 
tion a^, the effective angular velocity vector 
is still bounded by the surface of a circular 
cylinder of radius a^tV. 

The critical design loads in the positive 
and negative directions of the unit vector ^ are 
asymmetrical and different from the symmetrical 
design load in the direction of the unit vector 
j. The display format for the motion vector com- 
ponents given by equation (14) is shown in figure 
2(a), where the different critical design loads 


(a) Angular velocity. 



(b) Acceleration. 

Figure 2. — Display of mo- 
tion vector components. 
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are represented by different circles. However, the use of more than one cir- 
cle can be avoided, if desired, by using appropriate scaling factors for the 
different pointers. 

It should be noted that figure 2(a) displays the motion vector components 
that give rise to the maneuvering g forces. This has been done to facilitate 
pilot control of the motion components that produce these forces. Equation 
(10) gives the actual force components produced. The display format which 
gives the directions of these components can be obtained by rotating figure 
2(a) through a positive angle of 90°, as shown in figure 2(b). 

Implementation 

In practical applications, the g-load display can be mechanized to show 
the angular-velocity components as variable-length pointers emanating from the 
center of a circle whose radius is determined by the linear velocity of the 
aircraft. Inasmuch as this mode of presentation emphasizes the role of the 
angular velocity vector, while at the same time displaying the influence of 
the constraint imposed by the linear velocity, it is well suited to the appli- 
cation considered in this report. 

Figure 3 shows a diagram of 
the apparatus required to generate 
the constraint circle. Three body- 
mounted accelerometers provide 
three output signals whose values 
are equal to the three linear ac- 
celeration components a^, of 

the aircraft relative to bony axes. 
Also provided is an input signal 
which is Indicative of the maximum 
design load of the aircraft in 
which the apparatus is placed. The 
three acceleration components are 
then passed through three integra- 
tors which yield the three velocity 
components u, U, and w. The velo- 
Figure 3. — Apparatus for processing city components are squared and 

accelerometer data. summed to obtain . The remaining 

circuitry extracts the square root 

of V'^ , determines its reciprocal, and formulates the signal = a^jV which is 
indicative of the instantaneous critical angular velocity. This is the angu- 
lar velocity which would produce g loads sufficient to cause structural fail- 
ure at the given velocity. As seen in figure 3, the signal is then passed 
through a sine wave generator to an oscilloscope to produce the circle shown. 
Before proceeding to generate the instantaneous values of the angular-velocity 
components Q and R in accordance with equations (11) and (12), it is necessary 
to generate the angles a and g, and to measure the components p, q, v and the 
Euler angles 0 and <j> . 
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The velocity components generated 
by the apparatus shown in figure 3 may 
be used to determine the angles of at- 
tack and sideslip and their trigono- 
metric functions as shown in figure 4. 



The three body-mounted rate gyros 
shown in figure 5(a) are required to 
measure the instantaneous values pj 
and r of the vehicle's angular velo- 
city relative to body axes. The two 
position-gyros shown in figure 5(b) are 
used to measure the Eulerian angles 0 
and (f>. 

Figure 4. — Apparatus for generating 
This completes the apparatus needed angles of attack and sideslip and 
to measure and generate the information their trigonometric functions, 
required to formulate the components Q 
and i? in accordance with equations (11) 
and (12). The output from the pieces of 
apparatus described can be combined by 
the use of multipliers and summers, to 
yield the components of the inertial and 
gravitational forces acting on an air- 
craft during maneuvers. 




CONCLUSIONS 


(3-) Measuring body rates. 


A display has been designed to 
facilitate control of g loads during 
remotely piloted vehicle maneuvers and 
to provide an immediate indication of 
p-load constraint violations. It is 
shown that for an allowable accelera- 
tion a^, the angular velocity vector 
(i) is bounded by the surface of a cir- 
cular cylinder of radius a^lV, which has 
for its axis the linear velocity vector 
V. For a specified value of a^, the 
cylinder expands or contracts as the 
aircraft velocity varies, thus indicat- 
ing to the pilot the linear and angular 
velocities permitted for a given accel- 
eration constraint. In practical appli- 
cations it is adequate to generate a 
circular cross section of a cylinder. 

When this form of display is used to (b) Measuring Euler angles and gener- 

facilitate pilot control of g loads, ating their trigonometric functions, 

pointers move in a circle of radius 

a^/V and assume lengths indicative of Figure 5. — Apparatus for measuring 
the acceleration-producing components of rate and attitude. 
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angular velocity. An alternative mode of operation permits the pilot to 
select a display in which g loads are presented directly as variable length 
pointers emanating from the center of a circle whose radius is again deter- 
mined by the linear velocity of the vehicle. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, June 3, 1975 


10 



APPENDIX A 


CONDITIONS FOR THE PARALLELISM OF THE VECTORS V AND 87/ 

To determine the conditions for the parallelism of the vectors V and 
'dVI'dt, it is first necessary to show that 

97 ..fi A 

■^ = (m + VO + wk) 

Proof : 


Let 


u = (V ' i) 


therefore 


u = 


du _ d 
It 


(V • i) 


u = 




dV ^ - di 


dt 


i + V 


' It 


that is 


u = 


87 _ —1 A _ 


but 
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di. 

'dt 


9^ , - - A 

+ 0JX7, = (a X T, 


Substitution from equation (A3) in equation (A2) gives 


u = 


87 A_ _A 

• 1. + Cm^V) • i + V • ((Dx-/) 


This equation may be rewritten as follows: 
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87 — A _ A 

• 1, — (7 X 0)) *i+(7x[o/ ' i 


therefore 
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and 


therefore 


or 


V = 


SF A 

U * ^ 


w = 


'87 f 
u ■ *= 


AAA AA AA AA 

(iii + vj + wk) = ^ ’ (ii + 33 + kk) 
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(ui 


V3 + wk) = -Trr ' $ = — 


8t 


8t 


(A4) 


where $ is the idemf actor. 

The paralj^elism conditions must ensure that 9F/8t has no component per- 
pendicular to 7. Hence, it is sufficient to show that, subject to the estab- 
lished conditions 
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This equation may be rewritten in component form as follows: 
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(vw-wv)i + (wu~ULj )3 + (ui>-vu)k 


(A5) 


(A6) 


Moreover 


u = V cos a cos 3 
V = V sin 3 
w = V sin a cos 3 

Substitution of these values in equa_tion (A^) yields the following condi- 
tions for the parallelism of the vectors 7 and 97/9t 

\_(cosa sin Q cos Q)a - (sina)^l = 0 

(cos^ 3M = 0 

(sin a sin ^ cos ^)a + (cosa)'^ = 0 

Consequently, the condition of parallelism will be satisfied, provided 

that 
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